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Abstract 

We present the first application of phase field modeling to electrochemistry. A 
free energy functional that includes the electrostatic effect of charged particles leads 
to rich interactions between concentration, electrostatic potential, and phase sta- 
bility. The present model, explored for the equilibrium planar interface, properly 
predicts the charge separation associated with the equilibrium double layer at the 
electrochemical interface and its extent in the electrolyte as a function of electrolyte 
concentration, as well as the form expected of surface energy vs. potential ("elec- 
trocapillary" ) curves. 



1 Introduction 



The phase field technique has previously been applied to the time evolution 
of complex dendritic, eutectic, and peritectic solidification morphologies [1,2]. 
The present work was motivated by the mathematical analogy between the 
governing equations of solidification dynamics and electroplating dynamics. 
For example, the solid-liquid interface is analogous to the electrode-electrolyte 
interface. The various overpotentials of electrochemistry have analogies with 
the supercoolings of alloy solidification: diffusional (constitutional), curvature, 
and interface attachment. Dendrites can form during solidification and during 
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electroplating. Nonetheless, we find significant differences between the two 
systems. 

Electrochemistry has been chosen by the microelectronics industry as the de- 
position method for "copper damascene" interconnects on microchips because 
it is capable of filling trenches and vias [3,4] without "pinch-off" void formation 
problems associated with chemical or physical vapor deposition. Because plat- 
ing into such small features naturally involves large curvatures of the growing 
interface and large electric field gradients, a rigorous mathematical solution of 
the governing equations is warranted. 

Modeling the evolution of a sharp interface, particularly when the topology 
is changing during growth, is a challenging problem because the boundary 
conditions at that interface are dependent on the shape and location of the 
interface. The phase field method avoids these difficulties by introducing a 
new phase variable and an appropriate governing equation for that variable. 
The governing equations for the system can then be solved on a uniform grid, 
with the interface position as one of the results of the solution. A simple 
analytical electrochemical cell consists of a working electrode, an electrolyte, 
and a reference electrode. In the case of electroplating, there will also be a 
counter electrode. The phase field variable in this work describes the transition 
between the working electrode and the electrolyte. 

1.1 Length scales in electrochemistry 

There are three significant length scales in electrochemistry: 

(1) the thickness of the electrode-electrolyte interface, 

(2) the charge separation (of the capacitive double layer) and voltage decay 
distance, which are related to the concentrations of charged species and 
the dielectric constant, 

(3) and the long range concentration decay length due to diffusion and con- 
vection in the electrolyte (related to the ratio of diffusivity to interfacial 
velocity). 

In traditional modeling of the equilibrium electrochemical interface, the electrode- 
electrolyte interface is assumed to be sharp, as shown in Figure 1. The charge 
separation distance is considered diffuse in the electrolyte only, with any ex- 
cess charge on the electrode residing entirely at the metal surface. In the 
electrolyte, there may be some specific adsorption of ions on the interface, 
and a more diffuse distribution of ions further away. Some ions may be com- 
plexed with solvent molecules, giving them a larger effective size and limiting 
their approach to the electrode surface (we do not consider solvation in this 
work). The long range concentration decay length is absent for an equilibrium 
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Fig. 1. Schematic of the electrode-electrolyte interface. 

interface. This equilibrium regime is characteristic of models such as that of 
Gouy and Chapman. By requiring that the ions in solution must simultane- 
ously satisfy a Boltzmann distribution and Poisson's equation, a description 
of the charge distribution is obtained, with predictions for the surface energy, 
charge separation, and adsorption. 

In another limit of traditional modeling, primarily used for dynamics, the 
electrode-electrolyte interface is assumed sharp and the charge separation is a 
dipole layer. The long range concentration decay length is modeled by trans- 
port in a stagnant fluid boundary layer. In equilibrium, the concentration 
and electrostatic potential jump across the interface are related by the Nernst 
equation 

^ = ^° + ^m^ (1) 
zj- L> R 

where ip eq is the equilibrium electrostatic potential jump, ip° is the standard 
potential, R is the molar gas constant, T is Faraday's constant, T is the tem- 
perature, Co is the concentration of the electroactive species in the oxidized 
state, Cr is the concentration in the reduced state, and z is the number of 
electrons transferred in the redox reaction. A term —VmjK/zJ 7 is added to 
the right hand side of Eq. (1) to treat curved interfaces, where V m is the molar 
volume, 7 is the surface energy, and K is the curvature. In dynamic condi- 
tions, a chemical reaction rate description (such as the Butler- Volmer relation) 
is used to describe the relationship between potential jump and current. 

In the present phase field model, all three length scales are assumed to be im- 
portant (albeit different) and all are modeled together with the same equations 
in the metal and electrolyte. This simplifies treatment of complex geometries 
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and permits proper treatment of adsorption and its effect on interface kinet- 
ics. As usual for phase field models, though, the equations can be very stiff, 
requiring significant computational resources. 



1.2 Traditional double-layer theory (Gouy- Chapman) 



The best known sharp-interface model of the electrochemical interface is that 
of Gouy and Chapman [5,6]. In this model (and its variants), the electrode is 
not considered, the distribution of ions in the electrolyte is assumed to follow 
a Boltzmann distribution 



Cj = C°°exp 



RT 



(2) 



and to satisfy Poisson's equation 

P = Yl C i z i :F 

N 



= — e 



d 2 j) 
dx 2 



(3) 



Cj is the concentration of species j, C°° is the concentration of species j in the 
bulk electrolyte, Zj is the valence of species j, p is the charge density, e is the 
permittivity of the electrolyte, and the electrostatic potential field ip = ipoo in 
the far field. In a z:z electrolyte (cation and anion have equal and opposite 
valence), these equations reduce to 



dip /8RTC°°\ 1/2 . . 
— — = — ) sinn 

dx V e 



which has the solution 



tanh \zT (ip - -0oo) /4-RT] 
tanh [zF {fa - if);*,) /4RT] 



zT (if) - Ippo) 
2RT 



exp (-x/l^) . 



(4) 



(5) 



ipo is the electrostatic potential in the electrolyte at the electrode surface 
(x = 0) and l^ = [eRT '/ '(2C oc 'z 2 ? 2 )} 1 / 2 is the voltage/concentration decay 
length (the same as the Debye length). 

Using this model, the surface energy, surface charge, differential capacitance, 
adsorption and other interfacial parameters can be related to the voltage across 
the interface. The surface charge on the electrode a s is described by 



' (97 
dip 



(6) 



Mi 



where 7 is the surface energy, and subscript /^ denotes that the partial deriva- 
tive is taken at constant chemical potential, i.e., at constant activities (or 
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concentrations) of all species, as approximated on an inert mercury electrode. 
A plot of surface energy as a function of ipo is known as an "electrocapil- 
lary curve" and from Eq. (6) we see that its maximum corresponds to the 
so-called "potential of zero charge" , where there is no net surface charge on 
the electrode (and due to overall charge balance, none in the electrolyte). 



2 Phase Field Model 



2. 1 Components 



The phase field approach requires a set of equations that can describe the 
evolution of conditions at every point in the system. Thus we must identify 
the chemical components that describe the electrode and electrolyte simulta- 
neously. As our model problem, we consider a Cu metal electrode in contact 
with a CuSC>4 aqueous electrolyte. We choose our electrode to be a solid so- 
lution of Cu +2 and interstitial e~. We choose our electrolyte to be an aqueous 
solution of H 2 0, Cu +2 , and S07 2 . We assume that the partial molar volume 
of e~ is zero and that the partial molar volumes of Cu +2 , H 2 0, and SOj 2 are 
all the same, such that 

V m = J2V 3 X J =Vsj:X J (7) 

and 

EC; = ^> (8) 

s y s 

where V m is the molar volume, Xj = CjV m is the mole fraction of component j 
with partial molar volume Vj, N is the set of all components, and S is the set of 
substitutional components with partial molar volume Vs. This treatment of e~ 
as interstitials allows the motion of e~ to decouple from the other species. The 
simplification that all substitutional species have the same, constant partial 
molar volume allows us to consider diffusion separately from deformation. 



2.2 Free energy 



The free energy of this two phase system of iV charged species is given by 



F(<p,C u ...,C N ,il>)= I 

J V 



dV (9) 
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where F is the total Helmholtz free energy, fy is the Helmholtz free energy 
per unit volume, p is the phase field variable, and k v is the phase field gra- 
dient energy coefficient, tp — 1 denotes the electrode and tp — denotes the 
electrolyte. The first term in the integral of Eq. (9) describes the chemical 
energy, the second term describes the electrostatic energy, and the third term 
describes the interfacial energy of the phase field. 



2. 3 Equilibrium 



We compose a Lagrangian from Eq. (9) by requiring that mass of each species 
must be conserved, and that Poisson's equation (3) and the substitutional 
constraint Eq. (8) must be satisfied everywhere in the system. By taking the 
virtual variations of this Lagrangian, we find that the governing equations in 
equilibrium are 

Vj _ 

/I j — jy-Hs = constant, (fO) 



= — -/^Vy? — (VV0 , (11) 

and Poisson's equation (3). The electrochemical potential of species j is 

/'/ || • (12) 

The subscript s denotes the solvent, chosen arbitrarily from the set S. The 
complete derivation of these equations is given in Ref. [7]. 



2.4 Choice of thermodynamics 



To simplify our model as much as possible, we choose to describe the chemical 
Helmholtz free energy by ideal solution thermodynamics 

f v (p, C,-) = £ C j {(if + A^p (p) + RT In CjV m + W 3 g (p) } (13) 

N 

such that 

/x, = nf + Afi° jP (<p) + RT In CjV m + + W j9 (p) . (14) 

nf is the standard chemical potential of component j in the electrolyte, A/z° 
is the difference between the standard chemical potentials of component j 
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Fig. 2. Phase field double well and interpolation functions. 

in the electrode and in the electrolyte, and Wj is the barrier height of the 
phase-field double well for component j. Other choices are possible, but we 
select 

p (cp) = (6ip 2 - 15y? + 10) (15) 

for the phase field interpolation function p, which is a smoothed step function, 
and 

g{ V )= V 2 {l- V f (16) 

for the phase field double well function g. These functions are plotted in Fig- 
ure 2. 



2.5 Choice of parameters 



If we consider only the values of concentration and potential far from an 
equilibrium interface (<p — 0, <p — 1), we can make contact with the sharp 
interface picture. Analysis of Eqs. (10) and (11) yields the fact that the far field 
values of the electrochemical potentials Jij are equal. From Eq. (14) (evaluated 
at ip — and ip — 1), we see that, for all species, 

^-^ = ^m$C-^f. (17) 

These equations describe the relationship between the potential difference and 
the concentration of each component in the two phases. Normally in electro- 
chemistry, the electrode is considered to be pure metal and the equations are 
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only applied to the species involved in electron transfer. Evaluated in this 
situation, Eq. (17) is equivalent to the Nernst equation (1). 

From Eq. (17), we see that the standard chemical potential difference between 
the electrode and electrolyte A/i° can be broken into two terms 

^° = ^ ln ^M-^ M (18) 

Aip Rci = (ifj s — ifj L ^J and superscript Ref denotes the reference state for the 
standard potentials. The chosen value of the parameter Atp Rci can be shown 
to have a direct relationship with the surface energy of the electrochemical 
interface, the charge on the metal, and the deviation of the interface from 
the point of zero charge. As such, it is a materials parameter that is fixed for 
each metal/electrolyte system, even as the actual concentrations are changed. 
Aijj Rei is equivalent to the "Galvani" potential of the electrode-electrolyte 
interface, and is related to the difference between the work functions of the 
electrode and electrolyte. 

We choose as our reference state a metallic copper electrode in contact with a 
1 mol/L CuS0 4 aqueous electrolyte, as described in Table 1. In order to per- 
form equilibrium numerical simulations, some trace of H 2 and SOj 2 must 
be permitted in the electrode and some trace of e~ must be permitted in the 
electrolyte. We arbitrarily choose a mole fraction of 10" 6 for these trace compo- 
nents in the reference state. The remaining parameters are W e - = J/mol, all 
other Wj = 3.6 x 10 5 J/mol, k v = 3.6 x 10" 11 J/m, V s = 1.8 x 10" 5 m 3 /mol, 
and e = 6.95 x 10" 10 F/m. 

Fixing the reference concentrations, we can examine the consequences of vary- 
ing the potential difference around A-^ Ref . We simultaneously solve Eq. (17) 
for all four components, subject to charge neutrality and J2 N Xj = 1 in each 
phase. The resulting mole fractions of the electrode and electrolyte are shown 
in Figure 3. We can see that, over the entire range of Aip Ref , the electrode 
remains copper metal, with trace amounts of H 2 and SO J 2 . Over most of 

Table 1 

Values of reference concentrations. 





x^L,Ref 
X 3 


Cf' Rei /(mol/m 3 ) 




Cf Ref /(mol/m 3 ) 


Cu+ 2 


0.0180 


1000 


0.333 


55600 


SO4 2 


0.0180 


1000 


1.0 x 10" 6 


0.167 


e~ 


1.0 x 10" 6 


0.0556 


0.667 


111000 


H 2 


0.964 


53600 


9.99 x 10" 7 


0.167 
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Fig. 3. Mole fraction as a function of applied potential for the electrode-electrolyte 
equilibrium. Note that the slope of SOJ 2 in the electrode is twice that in the elec- 
trolyte. 



the displayed range, the electrolyte is predominantly water, with Cu +2 and 
SOj 2 obeying Nernstian behavior. 

The behavior of the trace amount of SOJ 2 in the electrode warrants discus- 
sion. If we assume the SOJ 2 concentration of the electrode is fixed, the Nernst 
relationship for SOj 2 indicates that its concentration in the electrolyte should 
move in the opposite direction from Cu +2 in the electrolyte, for a given change 
in ip s — ip L , because of the sign change on Zj. The requirement of bulk charge 
neutrality in the electrolyte prevents this from occurring, however. If the con- 
centration of Cu +2 in the electrolyte increases by an order of magnitude, charge 
neutrality requires that the concentration of SOj 2 in the electrolyte increases 
by the same amount (neglecting e _ ). In order for Eq. (17) to be satisfied for 
both Cu +2 and SOj 2 , i.e., for all species to be in electrochemical equilibrium, 
the concentration of SO4 2 in the electrode must then increase by two orders 
of magnitude. 

Another way to see the equilibrium relationship between these phases is on the 
quaternary phase diagram in Figure 4. The vertices of the diagram are pure 
Cu +2 , SOj 2 , H2O, and e _ . The plane shows the zero charge combinations 
and the tie lines between liquid and solid phases show the equilibria varying 
between Cu(M)-H 2 at negative electrode potentials and Cu(M)-CuS0 4 at 
positive electrode potentials. 
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Fig. 4. Phase diagram and charge-neutral plane for Cu(M)-CuS04-H 2 0. 
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Fig. 5. Profiles through the ID equilibrium interface for different values of Atp Rei . 
The shaded area denotes the interfacial region 0.1 < p < 0.9. 



2.6 Equilibrium profiles 



Eqs. (3), (10), and (11) are solved with finite difference techniques (to second 
order in space) on a uniform grid. The sharp, two-phase initial condition was 
allowed to relax to equilibrium. The numerical methods are discussed in de- 
tail in Ref. [7]. The ID interface between electrode and electrolyte is shown in 
Figures 5 and 6. The different curves are for fixed concentration and different 
values of Aijj . We see that, as one might expect, the phase field tp is not 
particularly sensitive to Aip Rci , but the charge distribution p is. All of the tp 
curves fall within 2% of a hyperbolic tangent profile, but the charge density 
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Fig. 6. Concentration profiles through the ID equilibrium interface for different 
values of Atp Rei . The shaded area denotes the interfacial region 0.1 < tp < 0.9. 



changes from a negatively charged electrode at Aifj = —0.12 V to a posi- 
tively charged electrode at A^ Ref = +0.12 V. At Aip Kcl = V, we see that 
there is essentially zero charge in the electrolyte, and the dipole charge in the 
electrode likewise sums to a very small charge (note that a dipole layer does 
not necessarily imply the presence of polar molecules [6]). 

The charge distributions in Figure 5 are a result of the ionic distributions 
in Figure 6. At Aip Rei = —0.12 V, there is a surplus of e - at the electrode 
interface. Cu +2 remains at essentially the electrode density to some distance 
into the electrolyte, indicating an adsorbed layer of Cu +2 (note the specific 
adsorption of S07 2 in each case). 

We examine the electrostatic potential profiles of Figure 5 in Figure 7. We 
fit each of the curves to = ip^ + (ipo — ?/>oo) exp (—x/l^) and find excellent 
agreement for the decay length with = 0.304 nm predicted by the Gouy- 
Chapman model. The deviation between fit and model at the left edge of 
the plot are because the Gouy-Chapman model assumes an abrupt electrode- 
electrolyte interface and the phase field model treats it as diffuse. 
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Fig. 7. Comparison of electrostatic potential profiles (light solid lines) with predic- 
tions of Gouy-Chapman model (heavy dashed lines). The shaded area denotes the 
interfacial region 0.1 < ip < 0.9. 

2. 7 Interface properties 

From a ID analysis of Eq. (9), the equilibrium surface energy is described by 



where the prime denotes the derivative with respect to position. The surface 
charge on the electrode is related to the surface energy by 



(this is identical with the sharp-interface prediction of Eq. (6)). 

We plot the surface energy, calculated with Eq. (19), for different values of 
Aip Rei in Figure 8. These curves are quite similar to the "electrocapillary" 
curves of surface energy measured on mercury electrodes [6] for different ap- 
plied potentials. In the electrocapillary experiment, the applied potential is 
varied at an inert electrode with a fixed Galvani potential. In our simulations, 
the Galvani potential is varied at a chemically reactive, but unbiased, inter- 
face. For each curve, there is a value of Aif) Rci for which the surface energy 
is a maximum. We see from Figure 9, which is calculated from either side of 
Eq. (20), that there is a maximum surface energy at the "potential of zero 
charge", just as is predicted by the Gouy-Chapman model. 




(19) 




(20) 
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Fig. 8. Surface energy as a function of A^ Ref . The lines are to guide the eye. 
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Fig. 9. Electrode surface charge as a function of Aip . The lines are to guide the 
eye. 

3 Conclusions 

• Phase field model of electrochemical interface exhibits double layer behavior. 

• Decay length of electrostatic potential are consistent with sharp-interface 
models. 

• Dependence of surface energy ( "electrocapillary curves" ) and surface charge 
on potential are consistent with sharp-interface models. 
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• Crucial presence of charged species in electrochemistry leads to rich inter- 
actions between concentration, electrostatic potential, and phase stability. 

• Dynamic treatments of this model arise naturally [7]. 
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